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Abstract 

We discuss a systematic method of analytically calculating the asymptotic form of quasi- 
normal frequencies of a four-dimensional Schwarzschild black hole by expanding around the 
zeroth-order approximation to the wave equation proposed by Motl and Neitzke. We obtain an 
explicit expression for the first-order correction and arbitrary spin. Our results are in agreement 
with the results from WKB and numerical analyses in the case of scalar and gravitational waves. 



^email: suphot@swu.ac.th 
^email: gsiopsis@utk.edu 



1 



The study of quasi-normal modes in asymptotically flat space-times has attracted a lot of 
attention recently [1-11] because the asymptotic form of quasi-normal frequencies was shown 
to be related to the Barbero-Immirzi parameter [12, 13] of Loop Quantum Gravity (see [14-18] 
and references therein). Quasi- normal modes form a discrete spectrum of complex frequencies 
whose imaginary part determines the decay rate of the fluctuation. The asymptotic form of high 
overtones is 



where Th is the Hawking temperature. This has been derived numerically [19-23] and subse- 
quently confirmed analytically [3,5]. The imaginary part is large making the numerical analysis 
cumbersome, but is easy to understand by noting that the spacing of the frequencies 27riTff 
coincides with the spacing of the poles of a thermal Green function on the Schwarzschild black 
hole background. The analytical value of the real part was first conjectured by Hod [24] based on 
the general form of the area spectrum of black holes proposed by Mukhanov and Bekenstein [25] 
[5 A — 4Gln/c, A; = 2, 3, . . ., where 5 A is the spacing of eigenvalues and G is Newton's constant 
in units such that h — c = 1) which is related to the number of microstates of the black hole 
through the entropy- area relation S = ^ A [26]. It has an intriguing value from the loop quan- 
tum gravity point of view suggesting that the gauge group should be 5*0 (3) rather than SU{2) 
(since we have k — ?> instead of A; = 2). Thus the study of quasi-normal modes may lead to a 
deeper understanding of black holes and quantum gravity. 

The analytic derivation of the asymptotic form of quasi-normal frequencies (1) by Motl and 
Ncitzke [5] offered a new surprise as it heavily relied on the black hole singularity. It is intriguing 
that the "unphysical" region beyond the horizon influences the behavior of physical quantities. 
Here we extend the results of ref. [5] by calculating the first-order correction to the asymptotic 
formula (1). We solve the wave equation perturbatively for arbitrary spin of the wave and 
obtain explicit expressions for the quasi-normal frequencies. Our results are in agreement with 
numerical results [21,6] in the case of gravitational and scalar waves as well as results from a 
WKB analysis [8] in the case of gravitational waves. 

The metric of a four-dimensional Schawrzschild black hole may be written as 



'n 



(2n + l)7ri + ln3 



(1) 



ds^ = -h{r) dt^ + 



dr"^ 



+ r'^dn'^ , h{r) = 1 



2GM 



(2) 



r 
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where G is Newton's constant and M is the mass of the black hole. The Hawking temperature 
is 

where Tq = 2GM is the radius of the horizon. A spin-j perturbation of frequency uj is governed 
by the radial equation 

l7/\ 1 iT7-/\,T, 



where y(r) is the "Regge- Wheeler" potential [27, 28] 

and L = + 1) is the angular momentum operator. We have j = 0, 1/2, 1,3/2, 2 for a scalar, 
Dirac, electromagnetic, gravitino and gravitational wave, respectively. We shall keep the discus- 
sion general. In fact, it will be necessary to avoid integer values of j throughout the discussion 
and only take the limit j — > integer at the end of the calculation. 

It is convenient to express eq. (4) in terms of the "tortoise coordinate" 

r* = r + ro In ( — - 1 ) (6) 
V^o J 

We obtain 

+ = (7) 

to be solved along the entire real axis. At both ends the potential vanishes (V^ — >• as r* — >• ±oo) 
so the solutions behave as ^ ~ giiwr* -p^^ quasi-normal modes we demand 

^ ~ e^^'"''* , ^ ±00 (8) 

assuming Reo; > 0. Let us express the quasi-normal mode as 

^ = e---/(r.) (9) 

Then f{r^) ~ 1 as r* -\-oo and /(r*) ~ g^"^^* as r* —00 (near the horizon). Following 
ref . [5] , we continue the coordinate r analytically into the complex plane and define the requisite 



3 



boundary condition at the horizon in terms of the monodromy of /(r*(r)) around the singular 
point r — ro. It is easily seen to be given by 

A^(ro) ^e-^'^'^"" (10) 

along a contour running counterclockwise. This will serve as the definition of the boundary 
condition at the horizon. The contour surrounding the singularity r — ro can be deformed in 
the complex r-plane so that it either lies beyond the horizon (Rer < Tq) or at infinity (r —>■ oo). 
Then the monodromy only gets a contribution from the segment lying beyond the horizon. It is 
convenient to change variables to 

z = a;(r* — inro) = a;(r + ro ln(l — r/ro)) (11) 

where we choose the branch that yields z — > as r — > 0. The potential can be written as a series 
in ^/z, 

which is also a formal expansion in 1 / ^/uJ. 

To compute the monodromy around the singularity r = tq, we shall deform the contour so 
that it gets mapped onto the real axis in the 2;-plane. Near the singularity 2; = 0, we have 
z —-^T^- We shall choose our contour in the complex r plane so that near r = 0, the positive 
real axis and the negative real axis in the z-plane are mapped onto 

argo; Stt argo; 
argr = 7r — , argr = ^ — (13) 

in the r-plane. These two segments form a 7r/2 angle (independent of argo;). To avoid the r = 
singularity, we shall go around an arc of angle 37r/2. This translates to an angle Stt around z — Q 
in the 2;-plane. 

By considering the black hole singularity (r = 0), we obtain two linearly independent solu- 
tions. They are of the form 

/±(r) = ri±^Z±(r) (14) 
where Z± are analytic functions of r. By going around an arc of angle of 37r/2, these change to 

/±(e3-/V) = e^'^(^±^>/V±(r) (15) 



This is an exact result which follows directly from the wave equation (7). To proceed further, 
we need to relate the behavior near the black hole singularity to the behavior at large values of 
r in the complex r-plane. We shall do this by solving the wave equation (7) perturbatively. This 
will enable us to write the wavefunction as a perturbation series in 1/ ^Juj. 
To zeroth order in 1/ ^Juj^ eq. (7) can be written as 

^-(^-) — 

where we replaced the potential V (eq. (5)) by its zeroth-order approximation (sec cq. (12)), as 
suggested in ref. [5]. The solutions arc given in terms of Bessel functions. They can be chosen 
to match the two exact solutions (14), respectively. We obtain 

/f(.)=e-v,(o) = e^/^J,,/,(.) (17) 
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whose behavior at infinity [z — > oo) is 

/f (^)~e-cos(^-7r(l±j)/4) (18) 

By our boundary condition, we need /(z) to approach a constant as ^ — > oo along the positive 
real axis in the 2;-plane. This is accomplished by adopting the linear combination 

/(o) ^ j(o) _ ^-^m fio) ^ e^-yji/W (^) (19) 

which approaches a constant at infinity, 

/(o) ^ _e-'^(i+i)V4 sin(7rj72) (20) 

on account of eq. (18). Going along the Sir arc around 2; = in the 2;-plane, we obtain 

/(o) (e3-^) = e^^^'+^^'/^ (/f {z) - e-'^^'!'' f^^ (z)) (21) 
whose behavior at infinity is found to be 

/(°) {z) - e-'^(i+^)^/^ sin(37rj/2) + e^^^-^)'l^ sin(27rj)e2*^ (22) 

using eq. (18) again. The monodromy (10) to zeroth order is deduced from eqs. (20) and (22) to 
be 
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which yields a discrete set of complex frequencies (quasi-normal modes) , 

^ = (2n + l)7ri + ln(l + 2 cos(7ri)) + o{l/^/n) (24) 

This reduces to the expression (1) originally obtained numerically in the case of even integer j 
(scalar and gravitational perturbations). It was derived analytically in refs. [3,5]. 
To go beyond the zeroth order, expand the wavef unction in 1/ -y/cJ, 

* = *(°) + ^^*W + o(l/u;) (25) 

The first-order correction obeys the equation 

^ + + l) = V^o6V^>^'' , SV{z) = ^ + ^ Viriz)) (26) 

as is easily deduced from eqs. (7), (12) and (16). The solutions of the first-order wave equation 
(26) may be written in terms of the solutions of the zeroth-order approximation (16) in a standard 
fashion, 

¥l\z) = C¥'\z) r ¥'^6V^^2^ - C¥'\z) r ^f6v¥2^ , C = (27) 
JQ Jo sm(7rj/2) 

where we are integrating along the positive real axis on the 2;-plane {z > 0). The large- 2; behavior 
is found to be 

poo 

¥l\z) r^c^±cos{z-7i{l+j)/4:)-c+±cos{z-7r{l-j)/4:) , c±±^C ¥l'^Sv¥l'^ (28) 

Jo 

To find the small- 2; behavior, we expand 

SViz) = + o{l/u) (29) 

where we used eqs. (12) and (26). It follows that 

^ z^^inc^ {z) + o( 1 /uj) (30) 
where G± are even analytic functions of z. For the desired behavior at infinity, define 



where ^ ~ o(l) and the dots represent terms of order higher than o{l/y/uj). Demanding ^ ~ e 
as z — > +00 fixes the parameter ^ to 

e = + , e+ = c++e'^^^/2 - c+_ , e- = c__e-'^^^/2 - c+_ (32) 

Then / = e'^^ approaches a constant, 

/(;,) ^ -e--(^+^)^/Sin(7rj72) |l - j (33) 

as z +00 along the real axis. Turning to values of z in the neighborhood of the black hole 
singularity (around z = 0), we observe that by going around a 37r arc, we have 

¥l\e^-'z) = e'^^''^^^'/^¥l\z) (34) 

Therefore, eqs. (15), (30) and (31) imply 

^{e^'^'z) = *(°)(e''^^z) - ieM^+j)i/^^L^ \'^^+\z) - e-^''^'/\¥l\z) - iC'^^!^\z))} (35) 

As 2; — > — 00 along the real axis, we deduce 

f(z) - e-'^(^+^')'/Sin(37rj72) |l - , ^ a\ + e^^^"^)'/^ sin(27rj) |l - , ^ b\ e^'' (36) 

[ y/-u;ro J t V=cJro J 

where 

e""^'/' (e+ + ^e- - e cot(37rj72)) (37) 

and S is a constant that can be easily calculated but is not needed in our discussion. The 
monodromy (10) to this order is 

M{ro) = -'^^^j^ |l + e-/2 _ + ecot(37ri/2))| (38) 

sm(7rj/2) 2^y-uJro J 

leading to the quasi-normal frequencies 

^ = (2n + l)7Tt + ln(l + 2 cos(7rj)) + ^^=^ (C- - i+ + i cot(37rj72)) + o(l/n) (39) 
J^H •\/n-|-l/2 

which includes a o{l/y/n) correction to the o(l) asymptotic formula (24). To obtain an explicit 
expression, we need the integral 

J{v,n)= dzz-^'^J,{z)JJz) = ^ )_ ^J, — (40) 
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On account of (17), we obtain 

c..^.^^(^+^) + ^-^'V(±,A±,/2) (41) 
12x/2 sin(7rj72) ^ ^ ' ' ■> ' > ^ ^ 

for the coefficients in the first-order wavefunction correction (28). After some algebra, we obtain 

from eq. (32) an exphcit expression for the combination of constants needed in the o{l/^/n) 

contribution to the quasi-normal frequencies (eq. (39)), 

^ , -/o^ n .^ 3^(^ + 1) + l-i^ sin(27rj) 

+ C cot(37rj/2) = (1 - ^) - 



where we used the identity 



24v^7r3/2 sin(37ri/2) 
xr2(l/4) r(l/4 + j/2) r(l/4 - j/2) (42) 



r(y)r(l -y) = (43) 

Eq. (42) has a well-defined finite limit as j approaches an integer. 
In the limit j 0"*" (scalar perturbation), eq. (39) reduces to 

^ = (2n + + ln3 + ^ r^(l/4) + o{l/n) (44) 

in agreement with numerical results [6]. 

For gravitational waves (j — > 2), we obtain 

g ^ (2„ + + ms . ^ '-^i^ r'(l/4) + own) (45) 

in agreement with the result from a WKB analysis [8] as well as numerical results [21]. Exphcit 
expressions for different spins may be readily deduced from the general expression (39) using 
(42). 

To summarize, we have derived analytically the first-order correction to the asymptotic form 
of quasi- normal frequencies (1) for four-dimensional Schwarzschild black holes. We solved the 
wave equation perturbatively based on the zeroth-order solution proposed by Motl and Neitzke [5] 
and obtained an explicit expression for arbitrary spin of the wave (eqs. (39) and (42)). Our results 
agreed with numerical results [21, 6] in the case of scalar and gravitational waves, as well as the 
analytical expression found by WKB analysis [8] in the case of gravitational waves. It would be 
interesting to extend our method to more general space-time backgrounds such as Kerr black 
holes [10] and higher dimensions [11]. 
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